Asymptotic quantum degeneracy at a 2D corner 
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In quantum mechanics, asymptotic degeneracy is often considered in the context of a particle 
in a symmetric double-well potential, and is the phenomenon whereby pairs of energy levels come 
together to form doubly degenerate levels in response to an increase in the separation, or depth of 
the two wells. Here we highlight a new kind of asymptotic degeneracy that can arise when a particle 
is bound to a surface formed by the intersection of two planes, when the intersection angle is > tt. 
To demonstrate this effect we consider the bound states of a charged particle subject its own 'image' 
potential in a highly polarizable wedge, as a function of the wedge opening angle, a. 

PACS numbers: 03.65.-w,03.65.Ge,03.65.Sq,41.20.Cv 
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Quantum mechanical degeneracy is commonly associ- 
ated with the invariance of the Hamiltonian under some 
symmetry group. Accidental degeneracies also exist, that 
apparently do not arise from any underlying symmetry 
group (although in such cases the group could just be 
difficult to find \V). Further to this, there are situa- 
tions when two energies come so close as to be prac- 
tically indistinguishable — asymptotic degeneracy. This 
kind of degeneracy can arise in situations when the prob- 
ability density is sharply peaked in two (or more) loca- 
tions in space, but where the probability of tunneling 
between these locations is negligible. Asymptotic degen- 
eracy attracted interest in the context of ID double-well 
potentials, where such systems were taken as a model for 
vacuum-vacuum tunneling in gauge theories [2]-[5] . More 
recently, there has been renewed interest in the double- 
well, as a system where aspects of the physics of Bose- 
Einstein condensates (BECs) may be observed [6 ; for 
BEC interferometry p" ; to generate mesoscopic entangle- 
ment between binary BECs [8]; and as a model system 
for a bosonic Josephson junction [9]. 

Here we consider a particle confined within a two- 
dimensional potential formed by the sharp intersection 
of two planar surfaces; explicitly the system to be inves- 
tigated consists of a charged particle that is bound — via 
its own 'image' potential, or Coulomb self-energy — to the 
surface of a highly polarizable material wedge, of open- 
ing angle, a (see figure IT]). We show that in such cases 
it is possible — through changing a continuously across 
its range, to 27r — to separate the wavefunction of the 
particle into two nearly independent parts, as in the case 
of a double-well, which display an associated degener- 
acy between pairs of energy levels. This is different from 
the usual process of splitting the wave-function, where 
the splitting is achieved through increasing the strength 
of the confining potential, or the spatial separation be- 
tween the wells. These results may also be relevant 



in the context of recent work investigating the proper- 
ties of electronic states bound to liquid helium surfaces 

(e.g. mm)' 



Suppose that a particle is subject to a ID potential, 
F(x), that is symmetric around a point, xq, and has 
two minima, one either side of xq (a double- well) . It 
has previously been observed that the eigenvalues of the 
Hamiltonian for such a system come in pairs that cor- 
respond to eigenfunctions that are either symmetric or 
anti-symmetric around xq. As the separation between 
the two potential minima is increased then these pairs of 
eigenvalues tend to coalesce at a rate that is exponential 
as a function of distance. Consequently, for large separa- 
tions, the system exhibits a degeneracy of symmetric and 
antisymmetric states O H] . The aforementioned splitting 
between pairs of energy levels is known to be given by, 
AE = {2h / m)ip^{xi^){dipi^/ dx)\xQ^ where ipi^ is an eigen- 
state of just one well, in the absence of the other [121 
§50]. 

We consider a generalization of this system to 2D polar 
co-ordinates (r, ^), where G [— 7r,7r], and the potential 
is symmetric around the line ^ = 0. The wave-function 
satisfies the additional boundary condition that it van- 
ishes within the region \6\ > a/2. If there is no other 
symmetry than V{r^O) = V(r, — ^), then there will be 
two types of states: symmetric, '0^+^(r, 0) = '0^+^(r, — ^), 
and antisymmetric, ilj^~\r.,6) = —ip^~\r^—6). Suppose 
that the wavefunction is confined by the potential so as to 
be sharply peaked along two channels, ~ ±Qf/2, that 
terminate where they meet, at r = 0. An approxima- 
tion to the states of this system can be given in terms of 
single well states, V^(^)(r, 0) c^ [V^L(r, 0) ± ^l)^{r, -0)] /V2 
(i.e. where ?/^l corresponds to a particle confined by the 
potential along the line 6 = -\-a/2). The difference in 
energy between the '^'^^^ states and the associated single 
well state, ?/^l, can then be estimated from the current 



operator, 



2m 



(1) 

Integrating ^ over the potential well associated with V^l, 
(i.e. the region, ^ > 0), and applying the divergence the- 
orem we obtain an expression for the energy difference, 



As an explicit demonstration of this effect we con- 
sider a charge bound within the 'image' potential of a 
highly polarizable material that forms a wedge of open- 
ing angle, a, when the particle is restricted from enter- 
ing the interior of the material (e.g. due to an electronic 
band-gap within the material that overlaps with the vac- 
uum [131 [14] — see figure [l]). 



V^^ 



m Jo 



V^lV^V^^^) - V^^^V^V^L dr = E^^^ - El, 



J6'=0 



where, V^ = {l/r)d/de. Along 6 = 0; V^(+)(r,0) = 
V2V^L(r,0); V^(-)(r,0) = 0; \/e^l^^^\r,0) = 0; and 
Ve^p''~\r,0) = V^V^V^L(r,0). Therefore the difference 
between the symmetric and antisymmetric energy levels, 
AE = E^~^ — E^^\ is related to the value of V^l in the 
following way (c.f. [3l|4]), 



AE: 



mi: 



V'L(r,0)(VeVL)(r,0)dr. (2) 



From Q observe that as '^l, or its gradient vanish along 
the line of symmetry — in this case the region around 
the corner where the two potential wells meet — the split- 
ting between symmetric and anti-symmetric energy lev- 
els goes to zero, and the E^~^^ and E^~'^ levels become 
degenerate. 



For a wedge composed of a material with permittivity 
62, and an opening of angle a, containing a charged par- 
ticle, and a second material of permittivity, ei, the elec- 
trostatic potential of the combined wedge+charge sys- 
tem has previously been given in [15]. This result is a 
generalization of the case where 62 ^ oo and 61 ^ 1, 
which was derived some time ago [16l [17] . We have ob- 
tained the image potential for arbitrary a through sub- 
tracting an integral expression for the potential of a point 
charge from the expressions given both in [15] and [17], 
remembering to divide the result by a factor of two (de- 
tails will be given in a subsequent publication). In gen- 
eral the resulting expressions are rather complicated, and 
for simplicity we consider the case when ei ~ 1, and 
r = (e2 — ei)/(e2 + ei) -^ 1, where the potential takes 
the comparatively compact form. 




ip{r,e) = 



167reor 



[kia)-f{e,a)] 



(3) 



FIG. 1: (Color online) Contour plots for the electrostatic 'im- 
age' potential given in (l3|, the dark red lines indicate the 
highest value of the potential (smallest absolute value), and 
the yellow, green, and blue lines indicate the potential's low- 
est values (largest absolute values). Figure (a) shows a wedge 
of opening angle, a = Stt/IO, and in (b), a = 37r/2. 



with, 



k{a) = 



f{0.a) = 



r d 
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(tt - a)(l - 7/^^/"+!)) - (tt + a){r] - 7/^/^) 



1 - 7^^/a 



I _ j^'27T/a 



1 + 7?27r/a ^ 2/^^/" cos(27r6>/a) 



This above result is valid for a three dimensional wedge 
that is infinitely extended along the z axis, with the as- 
sumed symmetry along z allowing us to neglect this co- 
ordinate throughout. The case of arbitrary F will be 
treated elsewhere. Figure [l] illustrates two contour plots 
of ([3|, showing the form of the potential when a < tt 
and a > TT. From the figure it can be observed that 
(|3| is an instance of the aforementioned case, where two 
potential wells extend along the lines d — ±a/2, and 
meet at r = 0. Note that expression ([3| has the correct 
limit in terms of localized images at the angles, a = tt/u 
{n integer): for example, at a = tt, k{a) vanishes iden- 
tically, and /(<9,a) = l/cos((9) (3.261(2) [13), so that 
the potential reduces to the known single image result of 
(f = — e/(167reo^), where y is the position on the vertical 
axis in figure [l] 
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FIG. 2: (Color online) The variation of the upper bounds on 
the energies obtained from (El, and (pi, as a function of a. 
The circles show the known analytic values of the first two 
surface states for a particle bound to a planar material with 
F^ 1. 



Neglecting spin and relativistic effects, the Hamilto- 
nian for a particle of charge e interacting with the elec- 
trostatic potential given in ([3| is. 
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r dr 
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dr 






+ Mr,^) (4) 



where (p{r^O) is given by (Is]). It is assumed that the 
eigenstates of (p| vanish along the surface, and inside the 
wedge {\0\ > a/2). We find upper bounds for the values 
of the energies of the first three states using the vari- 
ational method, with trial wavefunctions that are con- 
structed using the known hydrogenic limit at a = tt (e.g. 



ipo -^ Noye ^/^) as a guide. 



V^o(r,6>) = A^or^^^cos 



ttO 



-7o{0)r 



Mr.O) = A^ir^^sin 



. f27r0\ 

^ sin 



-fiiey 



\ a J 



Mr,e) = N2r^ 



a — r cos 



7tO 



cos — I e 
V a 



-ii^ey 



(5) 



where ^i{0) = Ui cos{pi7T0/a)-\-qi^ and the Ni (z = 0, 1, 2) 
are the normalization constants. The variational param- 
eters are given by, {mi^rii^Pi^qi}^ and are restricted so 
that they can take only positive values (it is also required 
that Pi < 1), and the constant, a, is chosen so that the 
orthogonality between tpo and ip2 is guaranteed. 



(mo + m2 + 2) 



h{{mi,ni,pi,qi}) 
l3{{mi,ni,pi,qi})' 



where, 



rV2 
Jo 



cos^{x)dx 



[jo{ax/7r) + 72(ax/7r)] 



mQ-\-m2+n ' 



At a = TT, we expect '02 to approach the first, rather 
than the second excited state. At this angle, the poten- 
tial possesses translational symmetry, V{x-\-5x) = V{x), 
rather than simply the refiection symmetry assumed for 
the other angles. Therefore the eigenstates should sat- 
isfy, iIj{x-\-Sx) -^ '0(x)e*^^'^^, which is not possible in the 
antisymmetric ('^i) case. The completeness of the eigen- 
states of the Hamiltonian guarantees the existence of the 
antisymmetric state for all other a. 

To carry out the variation of the energy we must cal- 
culate the integrals of ipiHipi and I'^^p over the interior 
of the wedge. It is possible to perform the integration 
with respect to the radial co-ordinate analytically and 
write the results in terms of gamma functions (not to be 
confused with the material contrast parameter, F, which 
is here set to 1) multiplying functions of 0. For example, 
the normalization of the state, t/^q, is given by. 



A^o 



22mo + %/Q. 



F(2mo + 2)72(^2 = ^o, ^2 = P2 = ^^2 = 0) ' 
and the associated variational ground state energy is, 
n^aN'^T{2mQ + 2) ^/^ cos2(x) 



E^ = {MH\^o) 



22mo+2^^ 

All ^70 



mo(mo + l/2) (mo + 1/2) 



I'r^' 



C\dx, 



where; A 



ml - {if- B = (2mo + 1)70 + 
(^)'[2tan(x)7^ + 7^'] + (meVS^eofi^)^; C = 7o' + 
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bility density close to the line 6> = 0, so that the system 
behaves as though the problem involved two independent 
half planes. Physically this can be understood as a re- 
duction in the average kinetic energy coming from the 

—ip{?i / 2mr'^){d'^ ip / dO'^) term in the energy density. 
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FIG. 3: (Color online) The probability distributions obtained 
from minimizing (V^il^lV^i). Figures (a)-(c) show IV^^I^ for 
i — 1,2,3 with a — 1.657 rad, while figures (d)-(f) show 
the corresponding distributions when a = 5.011 rad. Note 
that the scale of figures (c) and (f) is 3:1 compared to the 
remaining plots. 



vk) 1q '^ 9 — ^(^) ~ /{ax/n^a); and 70 and its deriva- 
tives are evaluated at ax/ir. Similar expressions hold 
for the quantities relating to the first and second ex- 
cited states. The values of these integrals were min- 
imized with respect to the parameters, {771^,71^,^^,(7^}, 
using a simplex routine ^19j. The same results were also 
obtained from a second program based on a different nu- 
merical library [20]. Figure [2] illustrates the resulting 
upper bounds on the energies, Eq, Ei^ and £^2, as a func- 
tion of the opening angle of the wedge, a. When a < tt, 
the three energies take distinct values that increase at 
differing rates, as a increases. However, as a approaches 
TT, and beyond, £^0 and Ei tend to the same value of 
—0.85 eV, remaining degenerate until the wedge becomes 
a half plane (a = 27r). We interpret this as an asymptotic 
degeneracy, as outlined above, and expect higher energy 
levels to pair up in a similar way. 

In conclusion we have shown that a corner can isolate 
two parts of the wave-function so that the system be- 
haves as two separate wells, with effectively degenerate 
symmetric-antisymmetric pairs of states. This is in 
contrast to the typical double-well system, where the 
analogous splitting of the states is directly related to 
the physical separation, or strength of the potential wells. 

To explain the origin of the degeneracy we use the dis- 
tributions associated with the energies of figure [2] as a 
guide. Figure [s] (a)-(c) shows these states in the cases 
when a < TT, and (d)-(f) when a > tt. Although there 
is no guarantee that these distributions will approach 
the exact probability distribution as the energy is min- 
imized, it is clear that the minimization of the energy 
beyond tt is achieved through a reduction in the proba- 
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